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Abstract: We show that the decay g — > t\bW~ is important and can even be 
dominant in the region of parameter space where it is kinematically allowed. We 
discuss phenomenological implications within the Minimal Supersymmetric Standard 
Model and models with broken R-parity. We consider the flavour diagonal case as 
well as a possible mixing between squarks of different generations. In the latter case 
also the decay g — ► t\c is potentially important. We show that the decay g —*t\b W~ 
is sensitive to the stop mixing angle. Furthermore we demonstrate that in scenarios 
with a higgsino-like LSP the gluino decays mainly into final states containing top 
quarks or a light stop if allowed by kinematics. 
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1. Introduction 

At the Tevatron as well as at the future Large Hadron Collider (LHC) the search for 
supersymmetric particles is among the main topics of their experimental programs. 
Here the strongly interacting supersymmetric partners of quarks and gluons, squarks 
and gluinos, are expected to have the largest cross sections. Their production as well 
as their decays have therefore been intensively studied in recent years 0, [|. 

In these studies it has been assumed that the gluino g decays either into qqi if 
kinematically allowed or into qqx®, q'qxf and gx® otherwise. Here q i} Xi > xf denote 
squarks, neutralinos and charginos, respectively. These decay modes have been used 
in the searches for gluinos (see e.g. |4j] and references therein). However, there exists 
also the possibility that the gluino decays via a three body decay into the lighter 
stop, namely: g — > tibW~. The necessary mass hierarchy m^m^ + m t > rrig > 

+ my/ + mb can be obtained e.g. in the minimal supergravity model as will be 
shown below. In the case that this mass hierarchy is realized in nature there are 
further gluino decays, violating flavour, into tic. To our knowledge this interesting 
possibilities have not been treated so far in the literature. 

In this paper we discuss the decays g — > t\ b W~ and t\c and possible signatures. 
Our framework is mainly the Minimal Supersymmetric Standard Model (MSSM) || 
with conserved R-parity. However, we also discuss possible implications of R-parity 
violation for the signatures of these two decay modes. The paper is organized as 
follows: in the next section we present our conventions as well as the formulas for 
the gluino decays. In Section |3| we present our numerical results for various scenarios, 
with and without R-parity. Finally in Section ^ we present our conclusions. 
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2. Conventions and the formulas for the widths 



The parameters relevant to our discussions are the soft susy breaking mass parame- 
ters for the squarks Mq { , M Vi1 M Di , the trilinear parameters A t>b , the gaugino mass 
parameters Mi ;2 , the gluino mass m g , the higgsino mass parameter /i and the ratio 
of the Higgs vacuum expectation values tan/3 = vijv\. Here i is a generation index. 

We give the formulas for complex parameters to be as general as possible although 
later in the numerical discussions we confine ourselves to real parameters to reduce 
the numbers of free parameters. It is well known that in the third generation left 
and right squarks mix due to the presence of the large Yukawa couplings. We give 
here for completeness the formulas for the mass matrix matrices as well as for the 
mass-eigenstates. Neglecting a possible generation mixing the mass matrices read 
as: 




M? LL M? RL 
Klr m Irr 



(q = b,t) 



(2.1) 



with 



M( LL = 

M Ill = 
MIr = 

t t RR 



m q,3 + m t + {\ ~ § sin 2 e w) cos 2/3 m| 
( M Irl)* = mt(A t -fJ.* cot/3) 
M u,3 + m t + | sin2 °w cos 2/3 m\ 
m q,3 + m l ~ i\ ~ | sin2 e w) cos 2(3 m| 

( M ^l) * = m b( A b ~ A** tan (3) 
m d,3 + m l - | sin2 ®w cos 2(3 m| 



(2.2) 

(2.3) 
(2.4) 
(2.5) 

(2.6) 
(2.7) 



The mass eigenstates qi are (51,^2) = {qL-,qR)Tl q with 



e tip « cos 9q sin 9q 



— sin 9 q e l(f ^ cos 9q 



(2-8) 



COS 9q 



\Klr\ 



y/\M$, LR \ 2 + (mi-Ml tLL y 
if q = arg(M (? - iLR ) . 



, sin6*g 



M Ill ~ K 



Mt LR ? + (ml - Ml LL f 



(2.9) 



The mass eigenvalues are 

m L = \ {^ M Ill + Ml RR ) =F ^(Ml LL -Ml RR y + 4\Ml LR \*) . (2.10) 
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Figure 1: Feynman graphs for the decay g — > t\ bW 



The decay g — > t\hW proceeds via top-quark exchange and sbottom exchange 
as shown in Fig. [l| The relevant part of the Lagrangian is given by: 

C = X a J r [c h lP L + c%P R ) ~g a ks + KsU (c\P L + 4Pr) g%s 

(2.11) 



+ C tb W;t^b + d l3 W; [ibm* - it*d% { j + /i.e. 

Here r, s are color indices whereas i,j denote the mass eigenstates. The are the 
Gell-Mann matrices with normalization *Yli a Tr(\ a ) 2 = 16. The couplings are given 
by: 

93 Rl 2 e-**3/2 Ai _ ( Q = bt ) 

(2.12) 



I Rn Vth , C+i 



V2 



Here 03 is the phase of the gluino mass parameter M3 and Vth is the (33) element of 
the CKM matrix. The partial width can be written as 

dT \C tb \ 2 



2567r 3 m?— — 
9 dsdt 



m 2 ) 2 + m 2 T 2 



^11 2 
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w 
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(2.13) 



with K(x,y,z) = y/(x - y - zf - Ayz, s = (p g - p^) 2 , t = (p § - p b ) 2 and u 
(p g —pw) 2 - The total width is obtained by integrating in the range 



(m w + m b ) 2 , s max = (m g - m^) 2 , 
m~ + m~ + m 2 ^ + ml — s (m| — m~ )(m 2 v — ml) 



' nuri.iimx 



T 



2 

K(s } m 2 g,m 2 h )K(s } m 2 w ,ml) 
2s 



+ 



2s 



(2.14) 



(2.15) 



We will also consider the possibility of g — > t\C. Using the approximate formulas 
for scalar top - scalar charm mixing as given in we get: 



with 



A; 



a. 



Y(g -> t x c) = — |e| m g \ 1 



21 



2 2 

— Id; 



16tt 2 



In 



M 2 \ V t * b V cb m 2 



(2.16) 



(2.17) 



m 2 ^ / 2m 2 v cos 2 /3 



(MJ )2 + M 2 Dt3 + M 2 Hi + A 2 b ) (2.18) 
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g 2 , f M 2 \ V.* h V cb m 2 h 
A R = Jt- In | 2 > t 4 2.19 

where is a high scale which we assume to be the Planck mass to get a maximal 
mixing. Mh 1 is the soft susy breaking Higgs mass term and Vtb and V c b are the 
respective elements of the CKM matrix. Assuming proper electroweak symmetry 
breaking one gets at tree level m 2 H = m 2 A0 sin 2 (3 — cos2/5m|/2 — |/i| 2 , where m^o 
is the mass of the pseudoscalar Higgs boson. The formulas for Al and A^ are the 
result of a single step integration of the corresponding RGEs assuming that the CKM 
matrix is the only source of flavour violation at the GUT scale. In this approximation 
the parameters Mq 2 , M|, 3 , M^, and Aj, can be evaluated at any scale because the 
induced error would be of higher orders. Therefore the expression should be treated 
as a rough estimate giving the order of magnitude for the mixing. For defmiteness 
we take the corresponding values of the parameters at the electroweak scale. In 
addition one should note that this approximation is an expansion in ra&/ (mw cos (3). 
Therefore one expects, that for small tan f3 the quality of this approximation is better 
than for large tan f3. 

In principle there is also the possibility that the gluino decays according to 
g — > t\u. However, this decay is suppressed by \V u b/V c b\ 2 — 1CT 2 in the approximation 
used above and will therefore be neglected in the following. 

Formulas for the decays g — > qqx®, g — > q'qXk aric ^ 9 ~^ 9x1 can De found in 0. 
Formulas for the cross section for gluino pair production as well as associated gluino- 
squark production including QCD corrections are given in |7J and for associated 
gluino-gaugino production including QCD corrections in p. 



3. Numerical Results 

In this section we present our numerical results for three different frameworks: (i) 
The MSSM without flavour violation, (ii) The MSSM with minimal flavour violation. 
In this case we assume that at an high energy scale the only source of flavour violation 
is given by the CKM matrix. RGE running induces in this case non- vanishing flavour 
violating couplings couplings between squarks, quarks and gluino. However, these 
couplings have to be small to respect the bounds on flavour changing neutral currents 



FCNCs 0. (iii) A scenario where R-parity is broken by bilinear terms fllOfl . This 



class of models resemble in many respects also the case of R-parity violation with 



trilinear terms |Tl[ violating lepton number as will be shown below. 

The parameter space relevant for our discussion is given by: rriq, + m t > 
rrig > + mw + va^. In the following we take the parameters freely without 
referring to a high scale scheme such as minimal supergravity (mSUGRA) or gauge 
mediated supersymmetry breaking (GMSB). However, we want to stress that this 
mass hierarchy can be obtained in mSUGRA. This can be seen by plugging the 
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following approximate solutions of the 1-loop RGEs (see e.g. [12 and references 
therein) in the inequalities above: m| — Mq + 6.2ML, m g ~ 2>.bM\/ 2 and m| = 



0.43M 2 + 4.55M 1 2 /2 + mf + 0.2M 1/2 A -M 1/2 y'2.25M 1 2 /2 + 1.13M 2 + 20.2m 2 /2 where 
Mx/2, M Q and v4 are the universal gaugino mass parameter, the universal scalar 
mass and the universal trilinear coupling at the GUT scale. 

In all numerical examples below we take nig = 500 GeV, mi = 660 GeV and 
for the first two generation squark mass parameters 600 GeV. Moreover, we fix = 
—865 GeV. The relevant Standard Model parameters used are: m t = 174.3 GeV, 
m b = 4.6 GeV, sin 2 6 W = 0.2315, a s {m z ) = 0.118 and a{m z ) = 1/127.9. The 
parameters m^, M D3 , tan/?, /i and M 2 define the various scenarios discussed below 
and are specified in Table |3TT| . We use the GUT inspired relation Mi = 5 tan 2 8wM 2 /3 
to reduce the number of free parameters. We will mainly discuss the dependence on 
the stop mixing angle as this is the most interesting parameter giving rise to the 
strongest dependence. The parameter A t can be calculated from mf and cos 6~ t and 
we have explicitly checked in all examples that \A t \ < 1 TeV avoiding problems with 
possible colour breaking minima. 



3.1 The MSSM without flavour violation 

In this section we discuss the case of the MSSM without any flavour violation. Before 
starting with numerical details, we want to note that in practice the T(g — ► tiWb) 
is very well approximated (within an error of 1% and below) by considering solely 
the top-quark exchange in Eq. ( |2.13j ) except for the parameter region where the 
lighter sbottom is nearly on-shell. The reason for this is that angular momentum 
conservation at the W-t-b vertex implies that either the t-b subsystem or one of the 
VU-boson squark subsystems form a P-wave. This in turn implies that the sbottom 
exchange is spin suppressed compared to the top-quark exchange. 



Table 1: Parameters for various scenarios. In all cases we have taken the squark mass 
parameters of the first two generations equal to 600 GeV, A b = —865 GeV, = 660 GeV 
and m-g = 500 GeV. We use the GUT relation for Mi = 5 tan 2 6 w M 2 /3. 



scenario 


m~ tl [GeV] 


M D3 [GeV] 


tan/3 


M 2 [GeV] 


H [GeV] 


m^o [TeV] 


1 


340 


580 


6 


150 


500 


1.4 


2 


380 


580 


6 


150 


500 


1.4 


3 


340 


580 


30 


450 


150 


0.86 


4 


340 


550 


6 


450 


150 


1.4 


5 


340 


550 


6 


1100 


600 


1.8 


6 


340 


580 


30 


1100 


600 


0.86 
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T(g -> [GeV] 




-1 -0.5 0.5 1 

cos 9i 



Figure 2: Partial decay width T{g — > i^W 6) as a function of cos#j for = 500 GeV, 
tan/3 = 6, mf 2 = 660 GeV, M Ds = 580 GeV, A b = -865 GeV, m~ tl = 340, 360, 380 and 
400 GeV. 



In Fig. we show the partial width T(g —>■ t\Wb) as a function of cos# 4 ~ for 
m% x = 340, 360, 380 and 400 GeV. One encounters a strong dependence on cos^- 
independent of the mass. This is due to the fact that the V^-boson couples only to 
left-handed fermions and correspondingly only to left sfermions. In case of cos &i = 
the lighter stop is a pure right state and the decay is only possible due to a "spin-flip" 
of the exchanged top quark. 

In Figs. |] and |] we show the gluino branching ratios as a function of cos Of 
for mi x = 340 and 380 GeV. In the first example the decay g — > t\Wb (full line) 
dominates for most of the range. We want to note that this is not a kinematical 
effect because mi + m b + m w ~ 2m t + m f o > m t + m h + m~± . In the second case 
the decay g — > t^Wb is about as important as g — > t&Xi" despite the fact that it is 
kinematically suppressed compared to the final states with the lighter chargino. In 
both cases the asymmetry in cos#£ is mainly due to interference effects in the width 
for g — > tbxf between stop and sbottom contributions. Together with information 
from stop production and decays [13| this asymmetry can be used to get information 
on the sign of cos Of. Note that this is asymmetry is not physical in the sense that 
it can be measured but is given by the structure of the theory implying that the 
experimental findings together with the consistency of the theory tells one which 
sign of cos^j is realized in nature. In addition we want to note that the branching 
ratio for the final state ttx? is maximal near cos = because in this case the right 
stop contributes most. We have checked that the behaviour shown hardly depends 
on tan f3. 

The situation is somewhat different in scenarios where the \y\ <C |M 12 as can be 
seen in Figs. B and R Independent of tan/3 the channels g — > t\Wb (full line) and 
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Figure 3: Gluino branching ratios for scenario 1 of Table |3.1| . In a) the lines correspond 
to g — > t\W~b + t\W + b (full line), g — > ttx? (dashed line) and 5 — > t&xr + ^Xi (dashed 
dotted line). In b) the lines correspond to g —> bb~X\ (dashed line), g — > 66x2 (dotted line), 
3 ^ Eg QVXi ( ful1 line ), 3 - > Eg 99X2 ( lon g short dashed line), and g -> Y, q>q > q^Xi+WXi 
(dashed dotted line), q and q' are summed over u,d,c,s. 




Figure 4: Gluino branching ratios for scenario 2. The parameters are specified in Table 3.1 



and in the text. In a) the lines correspond to g — > t{W b + i\W + b (full line), g — > ttxi 
(dashed line) and g — * tb~Xi + tbxf (dashed dotted line). In b) the lines correspond to 
9 -> bbxi (dashed line), ~g -> bb\% (dotted line), g -> J2 q QQXi (full line), <? -> Eg <?<?X2 
(long short dashed line), and g — > Egg' Q<fXi + QQ'Xi (dashed dotted line), g and are 
summed over it, d, c, s. 



(7 — > tbxi (dashed dotted) line show a strong dependence on cos % The reason is that 
the IV-boson couples the left states whereas the higgsino-like chargino couples mainly 
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Figure 5: Gluino branching ratios for scenario 3. The parameters are specified in Table |3TT| . 
In a) the lines correspond to g — > iiW~~b + t\W + b (full line), g — > ttXi (dashed line) and 
g — > tb~Xi + ibxx (dashed dotted line). In b) the lines correspond to g — > bbxi (dashed 
line), g — > bbx% (dotted line), g —>■ YlqQQXi (full line), g —> ^2 q qqX2 (long short dashed 
line), and g — > J2 q q ' QQ' Xi + Ql' Xi (dashed dotted line), q and q' are summed over u, d, c, s. 

a ) BR(<?) b) BR(<?) 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



cos 9i cos 9i 



Figure 6: Gluino branching ratios for scenario 4. The parameters are specified in Table |3~l| . 
In a) the lines correspond to g — > t{W~b + i\W + b (full line), g — > tixi (dashed line) and 
g — > tb~Xi + tbxf (dashed dotted line). In b) the lines correspond to g — > bb~Xi (dashed 
line), g — > 66^2 (dotted line), g — > J2qQ<jXi (full line), g — > ^2 q qqX2 (long short dashed 
line), and g — > Yl q q > QQ* Xi + Ql' Xi (dashed dotted line), q and q' are summed over u, d, c, s. 

to the right stop giving rise to the peak of g — > tbxf near cos^f = 0. The asymmetry 
in the sign of cos Of can be traced back to stop-chargino-bottom coupling which reads 
as — g R^Vji + YtR^Vji. Here the relative sign (in case of complex parameters relative 
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BR(g -> he + tic 




10" 



10" 



10" 



10-' 

-1 -0.5 0.5 

cos % 

Figure 7: Branching ratio g — > tic + iic as a function of cos#£ for scenarios 1 (dashed 
line), 2 (dashed dotted line), 3 (full line) and 4 (dotted line). The parameters are given in 
Table 13. II and in the text. 



phase) between and R\ 2 gives rise to the asymmetry. In both cases the final state 
tixi is less important because the lightest neutralino is to a large extent a down-type 
higgsino and the decay into the second lightest neutralino, which is mainly a up-type 
higgsino, is kinematically suppressed. This ordering is beside kinematics also the 
reason for BR(g — > bbx\) > Bi(g — > bbxX) m these scenarios. We want to stress here 
that for small tan (3 the gluino mainly decays into modes containing at least one top 
or a stop in higgsino-like scenarios as can be seen in Fig. ^|. 

The signature of g — > t{Wb clearly depends on the decay modes of the lighter 
stop. In the examples studied here t\ decays mainly into bxt and the chargino decays 
further mainly into xWq an d Xi^ u - Depending on the parameters chosen, other 2- 
body decay modes of t\ can become important []13) or higher order decay modes 
are important in case all the two-body tree-level decay modes are kinematically 
forbidden [H, |TJj]. If for example the decay t\ mode is t\ — > bWx\ one gets a final 
state g — > W + W~bbx\- The same final state can be obtained via the chain g — > 
ttxi — > W + W~bbxi- Similarly, one can show that for all stop final states one finds 
a gluino decay into a neutralino or a chargino that contains the same particles in 
the final state. However, in general the energy distribution as well as the angular 
distribution of the final state particles will be different, which is of course important 
for gluino searches and measurement of gluino properties. 

3.2 The MSSM with flavour violation 

Let us now turn to the case that case that the flavour violating coupling gluino - 
stop - charm-quark is non-zero. As mentioned in Section |2| we use the formulas 
given in |J to describe the mixing between top-squarks and the left scalar charm. 
In Fig. we display the branching ratios for the scenarios 1 - 4 of Table 3J.. We 
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Figure 8: Gluino branching ratio for scenarios 6 of Table |3.1| . The lines correspond to 
g — > t\W~b + t\W + b (full line) and g — > l\c + t\c (dashed line). 



have checked that the used values for the stop - scalar charm mixing are compatible 
with the bounds given in |J. One clearly sees that under the assumption of minimal 
flavour violation for small tan/3 the branching ratio is at most 1-2%. However, in 
case of large tan/3 this decay mode can be potentially large giving branching ratios 
up to 20%. In case that this decay is important the main consequence is a reduction 
in the multiplicity of the final state compared to the other gluino decays. 

As mentioned in Section |] the used approximation for the description of the 
scalar-charm - stop mixing is the result of a single step integration of the corre- 
sponding RGEs assuming the CKM matrix is the only source of flavour violation at 
the high scale. Therefore, the results obtained above should not be taken literally 
but as a demonstration of the expected order of magnitude under the assumptions 
above. Clearly, additional flavour violation in the squark sector at the high scale 
could enlarge the branching ratio for the decay g — *■ t\c. However, it cannot be much 
larger than 20% because otherwise the experimental bounds from low energy physics 
will be violated. Moreover, in case of additional flavour violation in the squark 
sector also the g — > t\U could be potentially large because there is not necessarily a 
suppression of the form \V u b/V c b\ 2 in such a case. 



3.3 R-parity violation 



In this section we are going to study scenarios where R-parity is broken. Here we 
focus on the model where one adds bilinear terms to the superpotential as well as to 
the soft SUSY breaking part JTOfl. This class of models can successfully explain the 



neutrino mass hierarchy as well as the neutrino mixing angles (see |15| and references 
therein). It also resembles in many respects the models with explicit lepton number 
violating trilinear couplings A and A' [IT]. This can easily be seen by a rotation of 
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the superfields where one transforms the bilinear terms in the superpotential into 
trilinear terms 1 leading to oc EifRj\i and A^- fe oc €iY^//j, |16|. Here Y E and Y D 
are the lepton and down-quark Yukawa couplings and q are the parameters of the 
bilinear terms violating lepton number in the superpotential. 

In what follows we study the scenarios discussed above as well as scenarios where 
the stop is the LSP and the gluino is the next heavier supersymmetric particle. The 
latter two scenarios can e.g. be obtained in GMSB scenarios [17| or in string scenarios 
. In the bilinear model the data from neutrino experiments imply relatively small 



R-parity violating couplings, e.g. |ej///| ~ O(10" 3 ) |TJj (for a discussion on gluino 
decays and production in the case of larger R-parity violating parameters see [|T^. PU| 
and references therein). This implies that branching ratios for the gluino decays in 
scenarios 1-4 are practically the same as in the R-parity conserving case because 
the R-parity violating decay modes have a branching ratio of at most O(10~ 7 ). How- 
ever, now the lightest neutralino will decay further giving rise to additional jets and 
leptons in the final state |2l], |22| compared to the R-parity conserving case (for recent 
discussions of neutralino decays with trilinear R-parity couplings see e.g. [p3fl ). 

In scenarios 5 and 6 of Table ^T] the stop is the LSP and the gluino is the next 
to lightest SUSY particle. Scenario 5 is a low tan/5 scenario whereas scenario 6 is 
a large tan f3 scenario. In the case of small tan f3 the three-body decay dominates 
practically with 100%. Moreover, the stop decays mainly into a lepton and a 6-quark 
p4}| . Thus, the signature is in this case 2 6-jets, a VT-boson and a charged lepton. 



In the case of large tan (3 also the decay into t±c becomes important as can be seen 
in Fig. [8[ The corresponding signature is in this case a 6-jet, a c-jet and a charged 
lepton. Note, that due to the Majoranna nature of the gluino the final state leptons 
can have same sign which clearly reduces the Standard Model background. 

Let us finally comment on the case of gluino LSP. It decays in these scenarios 
into the following final states: qqVi and q'ql^ 1 and gvi. For the scenarios discussed 
above the width of the gluino varies between O(eV) and O(keV). It turns out that 
in such a case final states containing top-quarks and/or bottom quarks dominate 
and that the branching ratios are sensitive to neutrino mixing angles similar to the 
case of neutralinos f22|. This, however, is beyond the scope of this paper and will be 



discussed in a dedicated paper [25 



4. Conclusions 

We have computed and studied the decays g — > t\W~b and g — > t±c. We have 
demonstrated that the branching ratio of g — > t\W~b is large and can even be 
dominant. This decay does not lead to new final states compared to the decays 
g — > ttXi or 9 ~ ¥ tbxf ■ However, it leads in general to different energy and angular 

1 This still leaves bilinear terms in the soft SUSY breaking part of the Lagrangian. 
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distributions of the final state particle compared to the decays into a chargino or 
a neutralino. This is clearly of importance for searches of gluinos at present and 
future colliders. In addition we have worked out possible signatures in models where 
R-parity is broken by lepton number violating terms. 
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